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Resumen 


El incremento urbano en zonas costeras ha conducido a la eutrofización, cuya acumulación de residuos orgánicos 
ha afectado también al estuario del golfo de Guayaquil; por ello, se determinó el estado trófico utilizando el índice 
Karydis, a una milla del margen costero de la provincia de El Oro, en tres zonas establecidas segun la afinidad 
hidrográfica. Además, se realizaron análisis multivariados para establecer la relación entre el microfitoplancton 
y las variables físico-químicas. Se evidenciaron diferencias significativas en niveles tróficos y estructura 
comunitaria por zonas. La zona A (canal de Jambeli), mostró mayor índice trófico, nitrato (3.8) y silicato (4.9), 
con menor salinidad, oxígeno disuelto y transparencia; la zona B (área de intercambio marino-costera), con 
mayor influencia oceánica, presentó valores tróficos medios (nitrato: 3.4 y silicato: 4.6) y la zona C (isla Puná), 
mostró mayor salinidad y los menores índices tróficos (3.2 de nitrato y 4.3 silicato). Las especies representativas 
fueron en la zona A Nitzschia longissima, Skeletonema costatum y Pseudonitzschia delicatissima; en la zona 
B: P. delicatissima, Guinardia striata y Leptocylindrus danicus, y la zona C Paralia sulcata, Thalassionema 
nitzschioides y S. costatum. Se determinó que el área de estudio es mesotrófica, no obstante, la zona A, presentó 
el mayor nivel trófico, debido probablemente a una mayor influencia de actividades humanas. 


Palabras clave. Análisis multivariado. Estuario. Eutrofización. Índice Karydis. Mesotrófico. 


Abstract 


Urban growth in coastal areas has led to eutrophication, which in turn has caused the accumulation of organic in 
the estuarine system of the gulf of Guayaquil. In this scenario trophic status was determined using the Karydis 
index in three zones established according to hydrographic affinity one mile from the coastal margin of the 
province of El Oro. Additionally, multivariate analysis were developed to establish the relationship between 
micro-phytoplankton and physical and chemical variables. Significant differences in trophic levels and community 
structure were found in each zone. Zone A (Jambelí channel), had the highest trophic index, nitrate (3.8) and 
silicate (4.9), and lowest salinity, dissolved oxygen, and transparency; zone B (coastal marine exchange area), 
the one with the greatest oceanic influence, had medium trophic levels (nitrate: 3.4 and silicate: 4.6); zone C 
(Puná island), had the highest salinity and lowest trophic indexes (3.2 nitrate and 4.3 silicate). The representative 
species in zone A were: Nitzschia longissima, Skeletonema costatum, and Pseudonitzschia delicatissima; zone 
B: P. delicatissima, Guinardia striata, and Leptocylindrus danicus; zone C Paralia sulcata, Thalassionema 
nitzschioides, and S. costatum. The study area was determined as mesotrophic. However, zone A presented the 
highest trophic level, probably due to a greater influence from human activities. 


Key words. Estuary. Eutrophication. Karydis index. Mesotrophic. Multivariate analysis. 
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Introduccion 


Elincremento de las poblaciones en las regiones aledafias 
a la zona costera ha sido motivo de preocupación por 
la problematica que implica el vertido de aguas de 
origen doméstico cuyo destino final es el mar (Borja 
y Dauer 2008, Flynn eft al. 2015). Uno de los efectos 
del permanente vertido de estas aguas residuales, es el 
incremento de nutrientes en la zona costera provocando 
la eutrofización. Los estudios sobre el estado trófico de 
las zonas costeras fueron iniciados por Nixon (1995), 
quien planteó la necesidad de buscar indicadores para 
medir la contaminación orgánica a través del nivel 
trófico (Cloern 2001). 


Dentro de los procesos de determinación del 
nivel trófico, se consideran las concentraciones de 
clorofila (Carlson 1977), nutrientes (Karydis 1992) 
o combinaciones de ambos (Vollenweider 1992). 
No obstante, los métodos más reconocidos utilizan 
indicadores biológicos como el microfitoplancton, 
pastos marinos y macroalgas (Livingston 2003), de los 
cuales, el primero es considerado el mejor indicador 
de los procesos costeros; además, su estructura 
comunitaria es determinante para conocer ambientes 
de distintos niveles tróficos (Tsirtsis y Karydis 1998, 
Garmendia et al. 2013). 


Los índices tróficos de Karydis y Trix, han sido muy 
utilizados en Europa y América, y han sido probados 
en ambientes tropicales como los estudios realizados 
en la zona costera del golfo de México (Medina-Gómez 
y Herrera-Silveira 2003, Herrera-Silveira et al. 2004 y 
Tapia-González et al. 2008) y en el Pacífico colombiano 
(Betancourt-Portela et al. 2011) que recomiendan 
aplicar índices para gestores gubernamentales en planes 
de manejo. 


La zona costera de Ecuador, presenta variabilidad 
estacional marcada con dos temporadas, la de lluvias 
(diciembre a mayo) y la época seca entre junio y 
noviembre (Twilley et al. 2001; Prado et al. 2015). En 
esta zona, se encuentra el golfo de Guayaquil (3?S, 
80°W) que es el más grande ecosistema estuarino en 
la costa del Pacifico Sudeste (12,000 km?). Comprende 
la línea de costa desde la parte sur de la provincia 
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del Guayas y toda la provincia de El Oro en Ecuador 
(Figura 1). Está formado por un estuario exterior que se 
inicia en el occidente de la isla Puná hasta el meridiano 
81° W, y un estuario interior que se extiende desde el 
oeste de la isla Puná hasta el borde costero conteniendo 
el Estero Salado, el río Guayas y el archipiélago de 
Jambelí (Cucalón 1996). 


En el golfo de Guayaquil se han realizado varios 
estudios sobre la composición y abundancia del 
microfitoplancton, entre los que destacan el de Jiménez 
(1983), quien identificó 189 especies de diatomeas 
y silicoflagelados; Cajas et al. (1998) determinaron 
que las aguas interiores del golfo de Guayaquil 
corresponden a un área de baja diversidad debido a 
la dominancia de Skeletonema costatum y Polymyxus 
coronalis; Gualancañay et al. (2003) señalaron que 
las diatomeas dominan tanto en época seca como 
en la lluviosa, y que los grupos de dinoflagelados, 
cianobacterias y flagelados se encuentran en menor 
abundancia; Prado et al. (2015) determinaron que en 
el estuario exterior del golfo de Guayaquil frente a la 
provincia de El Oro, dominan especies oceánicas como 
Rhizosolenia imbricata y Eucampia cornuta, mientras 
que en el estuario interior son más abundantes Paralia 
sulcata y S. costatum. 


El crecimiento sostenido en Ecuador durante los 
últimos 20 años, se basa principalmente en actividades 
de agricultura y de acuicultura, particularmente en la 
región sur del golfo de Guayaquil en la provincia de 
El Oro (INEC 2013). Al respecto Twilley et al. (2001) 
recomendaron la determinación de la calidad de agua 
de esta zona, pues la misma se habría deteriorado 
como consecuencia de la pérdida de grandes áreas 
de manglares, considerados filtros biológicos de 
importancia primordial (Turcios y Papenbrock 
2014), que luego fueron utilizadas para la instalación 
de granjas camaroneras. Sin embargo, no se han 
realizado estudios sobre calidad de agua del golfo de 
Guayaquil y zonas aledañas, considerando que los 
vertidos originados por las mencionadas actividades, 
tienen como destino final ese importante cuerpo de 
agua. A partir del programa de monitoreo hidrográfico 
desarrollado por el Instituto Nacional de Pesca, se 
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planteó la determinación del estado trófico y la 
respuesta de la estructura comunitaria, abundancia y 
diversidad del microfitoplancton a través de especies 
indicadoras de estos procesos, en la zona sur del golfo 
de Guayaquil frente a la provincia de El Oro. 


Material y métodos 


El área de estudio comprende el estuario interior del 
golfo de Guayaquil a una milla aproximada de la 
costa continental frente a la provincia de El Oro, en 
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donde se establecieron 17 estaciones de muestreo 
georreferenciadas con un GPS marca Garmin modelo 
GPSMAP 62sc (10 m de precisión). Con base en la 
afinidad de las condiciones hidrográficas y luego 
de realizar un análisis de conglomerados (Wards 
City block), se dividió el área en tres zonas: A: área 
adyacente al canal Jambeli (estaciones 1, 2, 3, 4, 17 y 
23), B: área de intercambio con aguas marino-costeras 
(estaciones 7, 8, 9, 11, 12, 13, 14 y 16) y C: zona 
interna al sureste de la Isla Puná (estaciones 15, 20 y 
21) (Figura 1). 
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Figura 1. Ubicación geográfica de las estaciones de muestreo en la zona costera de la 
provincia de El Oro entre junio y diciembre del 2012. 


Colecta y análisis de muestras 


Los muestreos se realizaron mensualmente desde 
junio hasta diciembre de 2012, correspondientes a 
la temporada seca de Ecuador, colectándose un total 
de 119 muestras, a nivel superficial entre las 08:00 
y 10:00 am. En cada sitio se realizaron mediciones 
de temperatura y salinidad con un Conductivity 
Temperature Depth (CTD) Seabird SBE-19 Plus V2 
(0.0001 de precisión). Simultáneamente se tomaron 
muestras de agua con botellas Niskin de 8 L de 


capacidad, para la determinación de nitrito, nitrato, 
fosfato y silicato, los cuales fueron analizados de 
acuerdo con la metodología de Parsons ef al. (1984). 
Además se efectuaron mediciones de transparencia 
mediante un disco Secchi de 30 cm de diámetro. 


De las mismas muestras obtenidas con botella Niskin 
se separaron alícuotas de 200 mL para el análisis 
cualitativo y cuantitativo de microfitoplancton, las 
cuales fueron preservadas con solución de lugol. Para 
la identificación de especies se utilizaron las claves de 
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Cupp (1943), Cleve (1951), Schiller (1971), Jiménez 
(1983), Tomas (1997) y Bérnard-Therriault et al. 
(1999). El recuento se realizó mediante el método 
Utermóhl referido en Hasle (1978), para lo cual se 
utilizaron cámaras de sedimentación de 50 mL y un 
microscopio invertido Leica DMi8. 


Análisis de datos 


Paraestablecer posibles diferencias zonales, estacionales 
y mensuales, se realizó el análisis de varianza no 
paramétrico de Kruskal-Wallis (Sheskin 2004), luego 
de comprobar el no cumplimiento de los supuestos del 
Anova como son normalidad (Kolmogorov-Smirnov) 
y homocedasticidad a través del índice de Bartlett 
(Zar 1996), expresándose los resultados en diagramas 
de cajas y bigotes de acuerdo con Boyer et al. (1997). 


La relación entre variables hidrográficas en cada zona 
se realizó mediante un Análisis de Componentes 
Principales (ACP) a partir de la matriz de correlación 
(Johnson y Wichern 1992, Legendre y Legendre 1998). 


Se determinó el estado trófico del área de estudio, 
utilizando los datos de nutrientes, con base en el índice 
de Karydis et al. (1983) y Karydis (1992) mediante la 
siguiente fórmula: 


IT= C/C-logXi +Log A, donde: 


IT: Índice trófico 

C: Sumatoria de la concentración de nutrientes por 
estación 

Xi: Concentración mensual promedio del nutriente 
por estación 

A: Nümero de estaciones 


La escala numérica que determina el nivel trófico 
es la siguiente: 


* <3: Oligotrófico 


* 3-5: Mesotrófico 
e 75: Eutrófico 
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Para establecer la relación entre las especies relevantes 
de cada zona con las variables hidrográficas, se realizó 
un Análisis de Correspondencia Corregido (DCA) y 
con base en la distancia del gradiente (<3) se efectuó 
un Análisis de Redundancia (ADR) (ter Braak y 
Verdonshot 1995, Muylaert et al. 2009). 


Las posibles diferencias en la estructura comunitaria del 
microfitoplancton entre zonas, se determinaron a través 
del Análisis Multidimensional no paramétrico (MDS- 
Anosim). Luego, se estableció la contribución de las 
especies relevantes del microfitoplancton por zonas 
mediante el SIMPER (Clarke y Warwick 2001). Para 
medir la respuesta de la comunidad ante variaciones 
entre diferentes tramos del área de estudio, se obtuvo 
la diversidad mediante los indices de Hill N1 y N2, y 
el de Shannon y Wienner (Tsirstsis y Karydis 1998, 
Chao et al. 2014). También se calcularon los índices 
de equitatividad y riqueza (Krebs 1999). 


Resultados 


Variables hidrográficas 


En la tabla 1 se detallan los valores promedio de las 
variables hidrográficas. 


Tabla 1. Promedios de las variables hidrográficas por zonas 
frente a la provincia de El Oro durante junio-diciembre de 
2012. Donde Prom = promedio y SE = Error estándar. 





Zona A B C 





Variable Prom. SE Prom. SE Prom. SE 
Temperatura°C 259 01 247 01 243 02 
Salinidad UPS 31,68 025 32,89 0,3 32,51 021 
Oxígeno Disuelto 6,02 024 727 026 69 023 
Nitrito umolL? 2,65 031 0,51 0,09 0,48 0,06 
Nitrato umolL? 7,18 0,83 2,51 033 947 08 
Fosfato umolL? 16 0,06 0,9 0,05 136 0,05 
Silicato umol L! 58,45 6,76 20,39 3,06 3494 6,87 
Transparencia(m) 0,9 0,5 255 025 12 0, 
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El comportamiento mensual de las variables 
hidrográficas fue bastante homogéneo durante los 6 
meses de estudio (p > 0,05). En cuanto a la variabilidad 
por zona todas mostraron diferencias significativas 
(p<0,05). La menor temperatura promedio (24,7 °C) se 
registró en la zona B, de mayor influencia costera y la 
más elevada (25,9°C) en la zona anexa al canal Jambeli 
(A) (Figura 2a). La menor salinidad promedio (31,68 
UPS) se detectó en zona A y la mayor (32,89 UPS) en la 
zona B (Figura 2b). La menor concentración promedio 
de oxígeno disuelto (6,02 mg/m?) se detectó en la zona 
A y la mayor (7,27 mg/m?) en la zona B (Figura 2c). La 
menor concentración promedio de nitrato (2,51 umol/L) 
se detectó en la zona B y la mayor (9,47 umol/L) en la 
zona C (Figura 2d). La menor concentración promedio 
de nitrito (0,48 mol/L) se detectó en la zona C y la 
mayor (2,65 umol/L) en la zona A (Figura 2e). La 
menor concentración de fosfato promedio (0,9 umol/L) 
se detectó en la zona B y la mayor (1,6 umol/L) en la 
zona A (Figura 2f). La menor concentración promedio 
de silicatos (20,39 umol/L) se detectó en la zona B y 
la mayor (58,45 umol/L) en la zona A (Figura 2g). 
La menor transparencia promedio (0.9 m) se detectó 
en la zona A y la mayor (2.6 m) en la zona B (Figura 
2h). En general la menor salinidad, concentración de 
oxígeno disuelto y transparencia se registraron en los 
alrededores del canal de Jambelí (zona A), así como 
elevadas concentraciones de nutrientes con excepción 
del nitrato que observó las mayores concentraciones en 
la zona C. En el área con mayor influencia oceánica (B), 
la concentración de nutrientes fue baja, no obstante, se 
reportó la mayor salinidad, concentración de oxígeno 
disuelto y transparencia. Mientras que la zona de la Isla 
Puná (C) alcanzó la menor temperatura y moderadas 
concentraciones de nutrientes y oxígeno disuelto. 


Al realizar el Análisis de Componentes Principales de 
la zona A, se observó el 65 % de varianza acumulada en 
los tres primeros ejes. En el primero, se correlacionaron 
en forma significativa la temperatura y los nutrientes 
(con excepción del nitrato) en relación inversa al 
oxígeno disuelto, salinidad y transparencia; mientras 
que en el componente 2 sólo se correlacionó con el 
nitrato (Figura 3a). 
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Figura 2. Variación espacial de las variables hidrográficas 
a) Temperatura; b) salinidad; c) oxígeno disuelto; d) nitrito; 
e) nitrato; f) fosfato; g) silicato, h) transparencia en la 
provincia de El Oro entre junio y diciembre del 2012. 
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El área norte de Jambeli mostró el 60 % de varianza 
acumulada en los tres primeros ejes, observándose 
correlación entre los nutrientes, los cuales estuvieron 
en relación inversa a salinidad y oxígeno disuelto. 
Mientras que en el componente 2, se detectó la relación 
entre temperatura y transparencia (Figura 3b). 


En la zona adyacente a Puná, la varianza acumulada 
fue de 71 % en los tres primeros componentes. En 
el primero, se correlacionaron temperatura con 
silicato y fosfato, en relación inversa a la salinidad y 
transparencia. Mientras que en el eje 2, se detectó la 
relación inversa entre oxígeno disuelto, nitrito y nitrato 
(Figura 3c). 
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Figura 3. Proyección ortogonal de los dos primeros 
componentes del ACP de las variables hidrográficas en a) 
Zona del canal de Jambelí, b) Zona de intercambio marino- 
costero y c) Zona adyacente a Isla Puná, entre junio y 
diciembre de 2012. 
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Niveles tróficos 


No se detectaron diferencias significativas mensuales 
en el nivel trófico en cada uno de los nutrientes, pero 
sí entre zonas (p < 0,05). Los valores promedio más 
elevados del índice trófico en nitrato y silicato se 
registraron en el área del canal Jambelí (zona A), de 
3,7 y 4,9 respectivamente. Los menores registros se 
observaron en el área con mayor influencia oceánica 
(zona B), con un valor promedio en nitrato y silicato 
de 3,1 y 4,2 respectivamente considerándose un nivel 
mesotrófico. Con excepción de nitrato en el que el 
mayor nivel trófico se registró en la zona adyacente 
a la Isla Puná (zona C) considerada mesotrófica; solo 
en silicato superó el nivel 5 (eutrófico) para el canal 
Jambelí (Tabla 2; Figura 4). 


No se detectaron diferencias significativas mensuales 
en el nivel trófico en cada uno de los nutrientes, pero 
sí entre zonas (p < 0,05). En el área del canal Jambelí 
(zona A) se registraron los valores promedio más 
elevados del índice trófico para todos los nutrientes. 
En la zona de intercambio marino costera (B), los 
valores obtenidos, fueron menores que en la anterior, 
en tanto que en los alrededores de la isla Puná (zona 
C), se observaron los menores registros (Tabla 2). 
Temporalmente, en junio y julio se detectaron elevados 
valores de índice trófico en silicato, los cuales fueron 
superiores a 5 en todas las zonas, sin embargo, a partir 
de agosto no superan valores de 4 (Figura 4). 


Tabla 2. Índice trófico de Karydis por nutrientes en tres 
zonas de la provincia de El Oro entre junio y diciembre del 
2012. Donde ITNO2: índice trófico para nitrito, ITNO3: 
índice trófico para nitrato; ITPO4: índice trófico para 
fosfato y ITSiO4: índice trófico para silicato. 














Zona ITNO2 ITNO: ITPO4 ITSiO4 
A 3,2 3,8 3,0 4,9 
B 2,5 3,4 2,9 4,6 
C 2,0 3,2 2,6 4,3 
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Figura 4. Variación mensual y por zonas del nivel trófico (Karydis) de a y b) nitrito, c y d) nitrato, e y f) fosfato, g y h) 
silicato, en la provincia de El Oro entre junio y diciembre del 2012. 


BIOTA COLOMBIANA 18 (2) - 2017 151 


152 


DOI: 10.21068/c2017.v18n02a09 


Microfitoplancton 


Se identificaron 183 especies distribuidas en 
145 diatomeas (en mayor proporcion del grupo 
Centricae), 33 dinoflagelados, tres silicoflagelados, tres 
cianobacterias y dos euglenophyceas (Anexo 1). Esta 
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estructura fue diferente para cada zona, registrándose 
el mayor número de especies en el canal de Jambelí 
(Zona A) y la menor en la aledaña a Isla Puná (Zona 
C) (Tabla 3). Las diatomeas, además de constituir el 
grupo con mayor número de especies, también fue el 
más abundante. 


Tabla 3. Distribución de los grupos funcionales del 
microfitoplancton por zona en la provincia de El Oro entre 


junio y diciembre del 2012. 























Grupo Zona A 
Diatomeas 133 
Dinoflagelados 33 
Silicoflagelados 3 
Cianobacterias 3 
Euglenales 3 
Total 105 


Zona B Zona C 
129 84 
26 17 
3 3 
3 1 
2 0 
163 105 





La abundancia promedio de microfitoplancton, no 
mostro diferencias significativas mensuales pero si 
entre zonas (p < 0,05) con el mayor registro en la zona 
B (Figura 5 a y b). La riqueza mostró comportamiento 
similar, con poca variabilidad mensual (p > 0,05) 





Abundancia log (x--1) 


Meses 





y diferencias significativas por zonas, con mayor 
promedio en la zona B (Figura 5 c y d). Mientras que 
la equitatividad no mostró variación significativa (p > 
0,05) mensual ni por zonas (Figura 5 e y f). 
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Figura 5. Variación mensual y por zonas de a y b) abundancia, c y d) riqueza, e y f) equitatividad 
del microfitoplancton en la provincia de El Oro entre junio y diciembre del 2012. 
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Cont. Figura 5. Variación mensual y por zonas de a y b) abundancia, c y d) riqueza, e y f) 
equitatividad del microfitoplancton en la provincia de El Oro entre junio y diciembre del 2012. 


Los registros de diversidad mostraron similar 
comportamiento, sin diferencias significativas (p > 
0,05) mensuales. No obstante, se detectaron diferencias 
entre zonas (p « 0,05) con los menores registros en la 
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zona del canal de Jambeli (A) y los mayores en la zona 
B que corresponde al área de intercambio con aguas 
marino-costeras (Figura 6). 
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Figura 6. Variación mensual y por zonas de: a y b) diversidad de Shannon, c y d) Indice de Hill 
N1, ey f) Índice de Hill N2, del microfitoplancton en las zonas A, B y C dela provincia El Oro 


entre junio y diciembre del 2012. 
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Cont. Figura 6. Variación mensual y por zonas de: a y b) diversidad de Shannon, c y d) Indice de 
Hill N1, ey f) Índice de Hill N2, del microfitoplancton en las zonas A, B y C de la provincia El 


Oro entre junio y diciembre del 2012. 


En la zona del canal de Jambelí, a través del Análisis 
de Redundancia (ADR) se observó en el componente 
1 la correlación entre Navicula transitrans, Paralia 
sulcata y Pleurosigma normanii con nitrato, así 
como Coscinodiscus oculus-iridis con oxígeno 
disuelto y salinidad. Todos ellos, en relación inversa 
a Pseudonitzschia delicatissima que presentó baja 
relación con nitrito. En el eje 2 Gyrodimium sp., se 
correlacionó con silicato y fosfato, en relación inversa 
con Skeletonema costatum, Chaetoceros decipiens 
y Nitzschia longissima asociada con transparencia y 
temperatura (Figura 7). 


En la zona de intercambio con aguas marino-costeras, 
en el componente 1, se correlacionó P. delicatissima 
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Figura 7. Proyección ortogonal de los dos primeros 
componentes del Análisis de Redundancia del 
microfitoplancton y variables hidrográficas en la zona A 
(Canal de Jambeli) frente a la provincia de El Oro, entre 
junio y diciembre del 2012. 


Troccoli-Ghinaglia y Prado-España 


con silicato, temperatura y fosfato, los cuales se 
presentaron en relación inversa a Thalassiosira subtilis, 
N. longissima y R. settigera. En el eje 2, Leptocylindrus 
danicus y S. costatum se correlacionaron con nitrito 
y transparencia, en relación inversa a Rhizosolenia 
imbricata correlacionada con salinidad y oxígeno 
disuelto (Figura 8). 
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Figura 8. Proyección ortogonal de los dos primeros 
componentes del Análisis de Redundancia del 
microfitoplancton y variables hidrográficas en la zona de 
intercambio con aguas marino-costeras frente a la provincia 
de El Oro, entre junio y diciembre de 2012. 


En la zona anexa a la Isla Puná (zona C), en el eje 1, se 
correlacionó Gyrodinium sp., con fosfato, y Chaetoceros 
decipiens, Coscinodiscus radiatus con silicato y 
temperatura, en relación inversa a Coscinodiscus sp., T. 
subtilis, Thalassionema frauenfeldianum asociados a la 
salinidad. En el componente 2 Pleurosigma normanii, 
Navicula transitrans y N. longissima asociadas en 
mayor grado a la transparencia del agua (Figura 9). 


La estructura comunitaria del microfitoplancton por 
zonas fue analizada a través del índice de similaridad 
(Anosim), denotando que a pesar que existen 
diferencias importantes (R: 0,253; p < 0,01) entre las 
tres zonas, también hay muchas especies comunes para 
toda el área de estudio (Figura 10). 
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Figura 9. Proyección ortogonal de los dos primeros 
componentes del Analisis de Redundancia del 
microfitoplancton y variables hidrográficas en la zona C 
(isla Puná) frente a la provincia de El Oro, entre junio y 
diciembre del 2012. 
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Figura 10. Diagrama MDS de la estructura comunitaria 
del microfitoplancton frente a la provincia de El Oro entre 
junio y diciembre del 2012. 
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Se registro variacion en la abundancia de las especies, 
es asi, que en la zona A se encontraron con mayor 
abundancia N. longissima, S. costatum, P. delicatissima; 
en la zona B fueron más abundantes las especies P. 
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delicatissima, Guinardia striata L. danicus, y N. 
longissima y finalmente en la zona C en mayor cantidad 
se registraron Paralia sulcata, T. nitzschioides S. 
costatum y N. longissima (Tabla 4). 


Tabla 4. Resultados del análisis de Simper del microfitoplancton en la zona 
costera frente a la provincia de El Oro durante la temporada seca de 2012. Donde 
Prom=promedio, % Contr=porcentaje de contribución, % Acum=porcentaje de 





acumulación. 
Zonas A B 
%Disimilaridad: 68 Prom Prom %Contr %Acum 





Pseudo-nitzschia cf. delicatissima 


197,54 317,84 9,85 9,85 





Nitzschia longissima 


330,61 287,63 6,83 16,68 





Guinardia striata 


62,61 167,22 5,48 22,16 





Leptocylindrus danicus 


30,18 171,13 5,45 27,61 





Skeletonema costatum 


171,34 125,71 5,43 33,03 





Zonas 
%Disimilaridad: 66 


A C 
Prom Prom %Contr %Acum 





Paralia sulcata 


71,1 306,06 9,54 9,54 





Nitzschia longissima 


330,61 207,28 8,1 17,64 





Pseudo-nitzschia cf. delicatissima 


197,54 18,9 6,66 24,3 





Skeletonema costatum 


171,34 173,96 5,17 29,46 





Thalassionema nitzschioides 


81,96 185,5 4,94 34,41 





Zonas 
%Disimilaridad: 66 


B C 
Prom Prom  96Contr %Acum 





Paralia sulcata 


33,44 306,06 9,06 9,06 





Pseudo-nitzschia cf. delicatissima 


317,84 18,9 8,17 17,23 





Nitzschia longissima 


287,63 207,28 6,33 23,56 





Guinardia striata 


167,22 52,61 4,93 28,49 





Skeletonema costatum 


125,71 173,96 4,83 33,32 





156 BIOTA COLOMBIANA 18 (2) - 2017 


Troccoli-Ghinaglia y Prado-España 


Discusión 

Al analizar el área de estudio por zonas, se determinó 
que los alrededores de la isla Jambelí (A), fue la zona 
más cálida, puesto que los promedios de temperatura 
fueron mayores en comparación con las otras zonas; 
esta variable y los nutrientes se correlacionaron 
significativamente, en especial el silicato que es un 
elemento conservativo (Day et al. 1989). Además se 
destaca la contribución de nitrito y fosfato en relación 
inversa al oxigeno disuelto, cuyas concentraciones 
fueron menores y que sirven de indicadores de 
actividad en los procesos biogeoquímicos en esa zona 
(Day et al. 1989; Van-Chu et al. 2014). Las variaciones 
en temperatura y salinidad, reflejan el intercambio entre 
masas de agua que resultan de la combinación entre 
los aportes fluviales y la zona costera adyacente en el 
área de estudio; estos intercambios son muy activos en 
ecosistemas con flujos de mareas significativos como 
ocurre en Ecuador (Fiedler y Talley 2006). 


El gran aporte de nutrientes en esta zona, probablemente 
es de origen continental, proviniendo de actividades 
que realizan las poblaciones ubicadas en áreas 
costeras de las islas. Gualancañay et al. (2003), 
indicaron que la zona sur del golfo de Guayaquil, 
recibe una mayor carga de nutrientes, debido a la 
gran actividad agrícola y acuícola que se realiza 
en este sector. La isla Jambelí actualmente cuenta 
con aproximadamente 1500 habitantes que realizan 
diferentes actividades relacionadas principalmente al 
turismo (e.g. restaurantes, viviendas) y acuacultura 
(INEC 2013), con el consecuente vertido de aguas 
residuales sin tratar, directamente al mar. 


La zona B, que es de intercambio entre las aguas 
marino-costeras, en la que se incluye los alrededores 
de la isla Santa Clara, registra la mayor oxigenación, 
salinidad y transparencia y menores concentraciones de 
nutrientes, notándose la influencia de aguas oceánicas 
que son más transparentes, así mismo, la correlación 
con el nitrato representa los procesos de oxidación de 
la zona costera adyacente, esto ha sido manifestado en 
varios estudios de la costa ecuatoriana, en los que se 
señala que la corriente de Humboldt, que se intensifica 
en la época seca, período en el que fue realizado este 
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estudio, proporciona una elevada carga de nitratos 
(Macías 1999), 


En la zona C, adyacente a la isla Puná, el análisis 
multivariado determinó una alta correlación entre las 
concentraciones medias de silicato y fosfato con la 
menor temperatura y además con valores medios de 
transparencia. Por su parte el nitrato registró los valores 
más elevados del área de estudio, debido a la influencia 
marcada de la corriente de Humboldt, que es fría y rica 
en nutrientes (Silva et al. 2009; Echevin et al. 2016). 
A pesar de ello, el nitrato no constituyó una variable 
relevante en el análisis de datos, por lo que su presencia 
no habría sido determinante en esta zona. 


Además de la presencia y disponibilidad de nutrientes, 
es fundamental considerar la luz como factor limitante 
en el crecimiento del microfitoplancton. En este 
sentido, aunque las tres áreas fueron turbias, la zona A 
presentó la menor transparencia (<1 m) durante todo 
el estudio; Pedersen et al. (2014) consideraron la baja 
transparencia como limitante para el desarrollo de la 
mayoría de las microalgas, sin embargo, las especies 
que dominaron como S. costatum, P. sulcata y T. 
nitzschioides son conocidas por proliferar en zonas 
con gradiente salino y en aquellas que están adaptadas 
a la intermitencia en la intensidad de la luz (Mitra et 
al. 2014). 


En cuanto a los niveles tróficos, los valores detectados 
fueron inferiores a tres de acuerdo a la escala de Karydis 
(1992), para el caso de nitrito, nitrato y fosfato, con 
excepción del silicato, que en meses específicos como 
junio y julio superaron valores de 5 como producto 
del aporte fluvial al golfo de Guayaquil, sin embargo, 
en promedio este nutriente no supera los niveles 
de eutrofización en la zona estudiada. Además, la 
estructura microfitoplanctónica estuvo dominada 
por diatomeas en aproximadamente el 90 %, lo que 
es un claro indicativo de la buena calidad del agua y 
justificaría la disponibilidad de silicato, el cual juega 
un rol fundamental en la presencia y abundancia de 
este grupo funcional (Egge y Aksness 1992, Wasmund 
et al. 2013). Al respecto, Lancelot y Muylaert 
(2011) señalaron la importancia de la relación entre 
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silicato y fosforo en los estuarios y zonas adyacentes 
que proporcionan condiciones favorables para el 
crecimiento de diatomeas. 


No seregistró abundancia de especies de cianobacterias 
que son indicadoras de ambientes eutrofizados; sin 
embargo, es necesario considerar que en la zona A, 
los valores del índice trófico son mayores que las otras 
dos localidades, lo que podría ser indicativo de que en 
esta área la producción de materia orgánica es mayor. 


Es importante señalar que en todas las estaciones, los 
muestreos se realizaron a una milla de la costa, en 
donde la influencia de las mareas es muy fuerte y por 
tanto determinante en la dilución de contaminantes 
orgánicos, sugiriéndose que el nivel trófico seria mayor 
en áreas cercanas al margen costero. En ese sentido 
Herrera-Silveira et al. (2004) manifiestan la necesidad 
de determinar el nivel trófico de localidades muy 
costeras, a distancias menores a una milla, para crear 
una base ambiental que permita establecer cambios 
espaciales y temporales del nivel trófico en sitios 
específicos y sobre todo en áreas de bahías y con poca 
circulación. 


Al asociar las especies de mayor frecuencia con las 
variables hidrográficas, se detectó la correlación de P 
sulcata con salinidad y silicato en las zonas de Puná 
y Jambelí, excepto en el área de intercambio con la 
zona costera anexa, donde destacaron especies de 
áreas marino-costeras para Ecuador como G. striata, 
P. seriata y L. danicus (Coello y Prado 1999). 


Al medir la respuesta de la comunidad ante variaciones 
ambientales, se observaron diferencias significativas 
entre zonas, registrándose una menor diversidad en la 
zona A, adyacente a la isla Jambelí. Las especies que 
dominan en áreas de baja diversidad, son resistentes 
a estresores (Muylaert ef al. 2009), y en el área de 
estudio fueron S. costatum, P. sulcata, T. nitzschioides 
y N. longissima que son organismos de alta tolerancia 
a cambios de salinidad, lo que las promueven sobre el 
resto de la comunidad a una alta abundancia, originando 
menor diversidad (Duarte ef al., 2006). 


Los Análisis de Redundancia mostraron correlación 
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significativa de todas las zonas en el primer componente; 
en la zona A, la relación entre altos valores de nutrientes 
y los menores registros de salinidad y oxígeno disuelto 
sugieren actividad en los procesos biogeoquímicos 
por los aportes fluviales (Howarth et al. 2011). Al 
respecto se ha señalado el papel del silicato al promover 
el crecimiento de diatomeas del grupo r-estrategas, 
es decir aquellas que tienen corto tiempo de vida, 
alta tasa de reproducción y adaptación a variaciones 
ambientales, como P. sulcata y S. costatum, forman 
parte del proceso de eutrofización y como consecuencia 
se pueden producir bajos valores en las concentraciones 
de oxígeno disuelto (Officer y Ryther 1980). En ese 
sentido, es importante señalar que los muestreos fueron 
diurnos; no obstante, en ambientes estuarinos, en horas 
de la madrugada se pueden alcanzar valores inferiores 
a 1 mg/L (Sawabini et al. 2015). 


Se registraron diferencias estructurales en la comunidad, 
con mayor heterogeneidad en la zona A, debido 
posiblemente a la interacción entre las aguas del 
golfo de Guayaquil y la zona costera adyacente que, 
durante la temporada seca, tiene influencia marcada de 
la corriente de Humboldt (Prado et a/. 2015). En ese 
sentido, se pudo establecer la contribución de especies 
relevantes, mostrando una clara diferenciación entre 
las diatomeas indicadoras de aguas costeras con las 
especies halotolerantes como P. sulcata, S. costatum, 
N. longissima y T. nitzschioides caracterizadas por 
tener früstulos gruesos. En tanto, que en la zona B, se 
presentaron P. delicatissima, G. striata y C. curvisetus 
que son especies con früstulos delgados y típicas de 
aguas marino costeras (Drira et al. 2014). Mientras que 
en la zona C (mas interna las especies representativas 
P. sulcata, T. nitzschioides y S. costatum son especies 
eurihalinas y características de ambientes estuarinos 
(Van Chu et al. 2014). 


Las especies S. costatum, P. sulcata, T. nitzschioides 
entre otras, han sido asociadas a procesos de 
eutrofización de ambientes estuarinos y costeros 
tropicales (Tas et al. 2009, Guo et al. 2014) por lo 
que se debe hacer seguimiento de la abundancia y 
distribución espacial y temporal para poder utilizarlas 
como indicadoras. 
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Otro aspecto que se debe considerar, es que la zona 
de estudio, tiene un margen costero con manglares, 
que puede jugar un papel importante en los ciclos 
biogeoquímicos de la zona, funcionando como filtros 
ante el flujo de nutrientes (Molnar et al. 2013, Herbeck 
et al. 2014) por lo que su presencia es fundamental para 
la salud del ecosistema. 


En general, se pudo establecer que de acuerdo a los 
valores de nivel trófico obtenido, la zona costera del 
golfo de Guayaquil, a una milla de distancia frente a 
la provincia de El Oro, es mesotrófica, no alcanzando 
niveles de eutrofización. Sin embargo, existen 
diferencias hidrográficas importantes entre zonas 
dentro del área estudiada, que se refleja en la estructura 
comunitaria del microfitoplancton y en los valores del 
indice trófico, el cual fue mayor en la zona A, lo que 
evidenciaría una mayor influencia de las actividades 
humanas en la isla Jambeli. 


La utilidad de aplicar índices es que si se prueba que 
funciona en ecosistemas costeros y estuarinos, puede 
usarse como herramienta para los gestores ambientales 
en los ámbitos costeros. Por ello se puede dar como 
ejemplo la utilización del índice Trix (multivariado) 
propuesto por Vollenweider (1992) y que se aplica en 
algunos países europeos. El índice Karydis, utilizado 
en esta zona del golfo de Guayaquil, constituye un 
aporte a la necesidad de establecer índices y probar su 
sensibilidad a nivel tropical. Al respecto, las variables 
para la detección de la calidad del agua son variadas 
y se pueden enumerar más de 300 a nivel global (Birk 
et al. 2012) y detectar los procesos de eutrofización 
debido al incremento de las actividades antropogénicas. 


Por otra parte, es necesario sefialar que, el presente 
estudio se realizó en la temporada seca por lo que 
se sugiere determinar los cambios estacionales en el 
estado trófico y de ser posible generar proyectos que 
incluyan la variabilidad interanual. 


Se recomienda realizar estudios en un área más 
cercana al margen continental en donde se involucre 
a los manglares en el área de estudio, puesto que la 
vegetación juega un papel importante en la filtración 
y asimilación de contaminantes procedentes de las 
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escorrentías en tierras altas, además de mejorar la 
calidad del agua, por lo cual es posible que la presencia 
de ellos haya sido determinante para que no existan 
niveles de eutrofización, especialmente en la isla 
Jambelí donde hay una gran influencia turística. 
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Anexo 1. Listado de especies del Microfitoplancton identificado frente a la Provincia EL Oro entre junio y diciembre del 


2012. 


DIATOMEAS PENNADAS 

Achnanthes sp. 

Amphora ovalis (Kutzing) Kutzing 

Amphora sp. 

Asterionellopsis glacialis (Castracane) F. E. Round 
Asterionellopsis kariana (Grunow) F. E. Round 
Bacillaria paxillifer (O.F. Müller) Hendey CC 
Climacosphenia moniligera Ehrenberg 

Cocconeis placentula Ehrenberg 

Cocconeis sp. 

Cylindrotheca closterium (Ehrenberg) 

Diploneis bombus Ehrenberg 

Diploneis bombus var. bombiformis (Cleve) Hustedt 
Diploneis sp. 

Diploneis splendida Cleve 

Entomoneis alata (Ehrenberg) Ehrenberg 


Fragilariopsis doliolus (Wallich) Medlin & Sims 
Gyrosigma balticum (Ehrenberg) Rabenhorst 
Gyrosigma fasciola (Ehrenberg) J. W. Griffith & Henfrey 
Lioloma delicatulum (Cupp) Hasle 

Lioloma pacificum (Cupp) Hasle 

Navicula brasiliensis Grunow 

Navicula cruciculoides C. Brockmann 

Navicula directa (W. Smith) Ralfs 

Navicula distans (W. Smith) Ralfs 

Navicula exigua (Gregory) Grunow 

Navicula lanceolata Ehrenberg 

Navicula sp. 

Navicula transitrans var. derasa (Grunow) Cleve 
Navicula transitrans var. derasa f. delicatula Heimdal 
Nitzschia acicularis (Kützing) W.Smith 

Nitzschia bicapitata Cleve 
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Nitzschia cf. sigma 

Nitzschia longissima (Brébisson) Ralfs 

Nitzschia navicularis (Brébisson ex Kutzing) Grunow 
Nitzschia obtusa W. Smith 

Nitzschia pacifica Cupp 

Nitzschia palea (Kutzing) W. Smith 

Nitzschia recta Hantzsch ex Rabenhorst 

Nitzschia sigma (Kutzing) W. Smith 

Nitzschia sp. 

Placoneis gastrum (Ehrenberg) Mereschkovsky 
Pleurosigma acuminatum (Kutzing) Grunow 
Pleurosigma angulatum (Queckett) W. Smith 
Pleurosigma angulatum var. genuinum (Queckett) W. Smith 
Pleurosigma delicatulum W. Smith 

Pleurosigma elongatum W. Smith 

Pleurosigma nicobaricum Grunow 

Pleurosigma normanii Ralfs 

Pleurosigma sp. 

Pseudo-nitzschia cf. delicatissima 

Pseudo-nitzschia pungens (Grunow ex Cleve) Hasle 
Pseudo-nitzschia seriata Cleve 

Thalassionema frauenfeldii (Grunow) Hallegraeff 
Thalassionema nitzschioides (Grunow) Mereschkowsky 
Ulnaria ulna (Nitzsch) P. Compere 


DIATOMEAS CENTRICAS 

Actinoptychus senarius (Ehrenberg) Ehrenberg 
Bacteriastrum delicatulum Cleve 
Bacteriastrum furcatum Shadbolt 
Bacteriastrum hyalinum Lauder 

Biddulphia alternans (J. W. Bailey) Van Heurck 
Cerataulina pelagica (Cleve) Hendey 
Chaetoceros affinis Lauder 

Chaetoceros borealis Bailey 

Chaetoceros brevis F. Schütt 

Chaetoceros coarctatus Lauder 

Chaetoceros compressus Lauder 

Chaetoceros curvisetus Cleve 

Chaetoceros debilis Cleve 

Chaetoceros decipiens Cleve 

Chaetoceros denticulatus Lauder 

Chaetoceros didymus Ehrenberg 

Chaetoceros didymus var. protuberans (H. S. Lauder) 
Chaetoceros eibenii Grunow in Van Heurck 
Chaetoceros laciniosus Schütt 

Chaetoceros laevis G. Leuduger-Fortmorel 
Chaetoceros lorenzianus Grunow 

Chaetoceros pendulus Karsten 

Chaetoceros peruvianus Brightwell 
Chaetoceros radicans Schütt 

Chaetoceros socialis H. S. Lauder 

Chaetoceros sp. 

Chaetoceros teres Cleve 

Chaetoceros tortissimus Gran 

Coscinodiscus asteromphalus Ehrenberg 
Coscinodiscus centralis Ehrenberg 
Coscinodiscus concinnus W. Smith 
Coscinodiscus curvatulus Grunow ex A. Schmidt 
Coscinodiscus granii Gough 

Coscinodiscus marginatus Ehrenberg 
Coscinodiscus oculus-iridis (Ehrenberg) Ehrenberg 
Coscinodiscus radiatus Ehrenberg 
Coscinodiscus sp. 
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Coscinodiscus wailesii Gran y Angst 
Cyclotella bodanica Eulenstein ex Grunow 
Cyclotella comta (Ehrenberg) Kiitzing 
Cyclotella meneghiniana Kützing 

Cyclotella sp. 

Cyclotella striata (Kutzing) Grunow 
Dactyliosolen fragilissimus (Bergon) Hasle 
Dactyliosolen phuketensis G. R. Hasle 
Detonula pumila (Castracane) Gran 

Ditylum brightwellii (T. West) Grunow 
Eucampia cornuta (Cleve) Grunow 
Eucampia zodiacus Ehrenberg 

Guinardia cylindrus (Cleve) Hasle 
Guinardia flaccida (Castracane) H. Peragallo 
Guinardia striata (Stolterfoth) Hasle y Syvertsen 
Hemiaulus hauckii Grunow ex Van Heurck 
Hemiaulus indicus Karsten 

Hemiaulus membranaceus Cleve 

Hemiaulus sinensis Greville 

Lauderia annulata Cleve 

Leptocylindrus danicus Cleve 
Leptocylindrus mediterraneus (H. Peragallo) Hasle 
Leptocylindrus minimus Gran 

Lithodesmium undulatum Ehrenberg 
Melosira moniliformis (O. F. Muller) C. Agardh 
Odontella aurita (Lyngbye) C. Agardh 
Odontella dubia (Brightwell) Cleve 
Odontella longicruris (Greville) Hoban 
Odontella mobiliensis (Bailey) Grunow 
Odontella sinensis (Greville) Grunow 
Paralia sulcata (Ehrenberg) Cleve 
Planktoniella sol (Wallich) Schutt 

Proboscia alata (Brightwell) Sundstróm 
Pseudosolenia calcar-avis (Schutze) 
Rhizosolenia bergonii H. Peragallo 
Rhizosolenia hebetata J. W. Bailey 
Rhizosolenia hyalina Ostenfeld 

Rhizosolenia imbricata Brightwell 
Rhizosolenia ostenfeldii B. G. Sundstrom 
Rhizosolenia robusta G.Norman ex Ralfs 
Rhizosolenia setigera Brightwell 
Rhizosolenia styliformis Brightwell 
Skeletonema costatum (Greville) Cleve 
Stephanopyxis palmeriana (Greville) Grunow 
Stephanopyxis turris (Greville) Ralfs 
Thalassiosira angulata (W. Gregory) Hasle 
Thalassiosira eccentrica (Ehrenberg) Cleve 
Thalassiosira lineata Jousé 

Thalassiosira sp. 

Thalassiosira subtilis (Ostenfeld) Gran 


CIANOBACTERIAS 
Anabaena sp. 
Leptolyngbya sp. 
Oscillatoria sp. 


DINOFLAGELADOS 

Ceratium dens Ostenfeld & Schmidt 
Ceratium falcatiforme E. G. Jorgensen 
Ceratium furca (Ehrenberg) 

Ceratium fusus (Ehrenberg) Dujardin 
Ceratium lineatum (Ehrenberg) Cleve 
Ceratium massiliense (Gourret) 
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Ceratium teres Kofoid 

Corythodinium sp. 

Dinophysis acuminata Claparede 
Dinophysis caudata Saville-Kent 
Diplopelta asymmetrica (Mangin) Balech 
Gonyaulax polyedra Stein 

Gonyaulax polygramma Stein 
Gonyaulax sp. 

Gymnodinium catenatum Graham 
Gyrodinium sp. 

Gyrodinium spirale Kofoid & Swezy 
Noctiluca scintillans Kofoid & Swezy 
Ornithocercus sp. 

Ornithocercus steinii Schutt 

Oxyphysis oxytoxoides Kofoid 
Oxytoxum sp. 

Prorocentrum gracile Schutt 
Prorocentrum micans Ehrenberg 
Protoperidinium conicum (Gran) Balech 
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Protoperidinium depressum (Bailey) Balech 
Protoperidinium elegans (Cleve) Balech 
Protoperidinium latispinum (Mangin) Balech 
Protoperidinium pentagonum (Gran) Balech 
Protoperidinium sp. 

Protoperidinium spirale (Gaarder) Balech 


DINOFLAGELADOS 
Pyrocystis noctiluca Murray ex Haeckel 
Pyrophacus steinii (Schiller) Wall & Dale 


EUGLENOPHYTA 
Phacus sp. 
Phacus suecicus Lemmermann 


SILICOFLAGELADOS 

Dictyocha fibula Ehrenberg 

Ebria antigua var. simplex Schultze 
Octactis octonaria (Ehrenberg) Hovasse 
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